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Abstract 

A generalized Lyapunov method is outlined which 
predicts global stability of a broad class of dissipa- 
tive dynamical systems. The method is applied to 
the complex Lorenz model and to the Navier-Stokes 
equations. In both cases one finds compact domains 
in phase space which contain the u sets of all tra- 
jectories, in particular the fixed points, limit cycles, 
and strange attractors. 

1 Introduction 

In the theory of ordinary differential equations, the 
method of Lyapunov function in general serves to 
examine the stability of a fixed point and its do- 
main of attraction, for an overview see e.g. jij]. 
The method can be naturally extended to the case 
where, instead of a fixed point, one is interested in 
the stability of a compact domain which has finite 
measure in phase space and is invariant under the 
flow. 

In this contribution we essentially illustrate the 
power of generalized Lyapunov functions, which are 
not discussed in the standard textbooks on dynam- 
ical systems theory. In particular we are interested 
in finding globally attractive domains for a certain 
class of nonlinear models. The corresponding sys- 
tems turn out to be globally stable in the sense that 
no trajectory which starts within a certain domain 
can leave it and the trajectories which start outside 
of the domain will end in it after sufficiently large 
times. The method presented cannot give details 
on the nature of the attractors contained within an 
attractive domain. Furthermore, we will only par- 



tially succeed to determine minimal attractive do- 
mains. This, on the other hand, opens the chance 
of finding attractive domains in an analytical way. 
As a matter of fact, in his famous paper, Lorenz 
showed that (nonminimal) attractive domains 
can be found in an elementary way by linear meth- 
ods provided the nonlinearities of the dissipative 
dynamical system are quadratic only and do not 
contribute to the overall energy balance. 

In the next section the method of generalized 
Lyapunov functions will be introduced together 
with a class of dynamical systems as proposed by 
Lorenz Q, which allow for quadratic Lyapunov 
functions. As a first example, the method is ap- 
plied to the real Lorenz model according to ||]. We 
discuss then in the third section a more detailed ap- 
plication to the five-dimensional or complex Lorenz 
model H . As compared with recent work || where 
the attractive domain has been successfully mini- 
mized to some extent for parameter values relevant 
in infra-red laser physics, further new results @ are 
presented here. In the last section the generalized 
Lyapunov method serves to prove the boundedness 
of the velocity field of the incompressible Navier- 
Stokes equations. This is a known result and was 
shown in different ways elsewhere, for the case of 
periodic boundary conditions see e.g. section 5.3 in 



2 Generalized Lyapunov func- 
tions 

Let us consider an autonomous dynamical system 

dx 
~dt 



R r 
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where the vector field / is supposed to be suffi- 
ciently smooth. Be L(x) a positively definite, suffi- 
ciently smooth, scalar function with L(0) = and 
L(x) > for x ^ 0. Furthermore, be G a domain 
which contains the point x = 0, and G its com- 
plement. Then we call L a generalized Lyapunov 
function, if the following properties hold for a com- 
pact domain G 

L(x) > for x / and L(0) = (2) 



dL d 

_:=-L(x(t,x„)) 



1=1 ° Xl 



for x G G, 
grad x L(x) ^ for x^O, 
and Lipschitz continuous. 



(3) 
(4) 



The manifolds L(x) = C = constant are closed 
hyper surf aces, which surround the point x = and 
foliate the phase space. To see this, one starts with 
a sufficiently small constant C which because of 
L(x) > is connected with an ellipsoidal surface. 
Then one constructs the one-dimensional curves 
normal to the surfaces L(x) = C. If A is a suitable 
curve parameter, the curves x(X) with x : R — > R n 
can be defined through the property that the curve 
tangents are parallel to the surface normal in every 
point. We can thus consider the curves as trajecto- 
ries of the following dynamical system 

dxi _ dL(x{\)) 
dX dxj ' 



l,2,...n. 



(5) 



Since this system is autonomous and because of (4) , 
we obtain unique curves x(X,xq) which do not in- 
tersect or touch each other; xo € R n is an arbitrary 
initial point which can be chosen, for instance, on 
an ellipsoid close to the origin. In this way we have 
constructed a 1-1 map between the points of an ar- 
bitrary surface L = C and the points of an ellip- 
soid close to the origin x = 0. In other words, the 
surfaces L = C are homeomorph to an ellipsoid sur- 
rounding x = 0. Moreover, any surface with a given 
constant Co separates the phase space into an in- 
ner part which contains x = and an outer part 
foliated by the surfaces C > Cq. 

At the boundary of G we have points where 
dL/dt = 0. We remark that this set contains 



the critical points of the system where f(x) = 0, 
because dL/dt = Y^jidL/dxi = 0. The surface 
L = C* which both bounds a domain containing 
the points with dL/dt = and is minimal with 
respect to C, is called the critical one. Now we 
are ready to draw conclusions for the trajectories of 
the system (1). If xq is the initial point of a tra- 
jectory with the property L(xq) = Cq > C* , then 
because of (3) we have dL{xo)/dt < 0. Therefore 
the trajectory wanders towards inner points with 
smaller C until the critical surface with L(x) = C* 
is reached. This tells that the critical surface is 
attractive from the outside. Simultaneously there 
can be no escape of a trajectory which starts inside 
the domain bounded by the critical surface. For 
illustration see Fig.l. The existence of a general- 
ized Lyapunov function guarantees therefore that 
the trajectories of the dynamical system asymptot- 
ically are confined to the domain G* bounded by 
the critical surface. Clearly, if we find several Lya- 
punov functions with attractive domains G*, G\, .., 
then the intersection G* = HiG* contains the mini- 
mal attractive domain. 




Figure 1: Illustration of a generalized Lyapunov 
function. The dashed curve shows the surface 
dL/dt = 0, which is tangent to the critical sur- 
face L = C* (dot-dashed). The latter confines a 
domain of attraction. The solid curves refer to sur- 
faces L = C with C >C* and dL/dt < 0. 

As a rather general example Lorenz considered 
the following dynamical system, for i = 1, . . . n, 

dxj 



dt 



^ ] Q>ijk%j%k ^ ] bijXj + Cj (6) 



3,k 



2 



with 

^2 a ijk XiXjX k = and ^x i b ij x j >0. (7) 

i,j,k i,j 

He proposed the Lyapunov function 

L{x) = ^{x 2 + x 2 + ...x 2 n ) (8) 
which gives rise to 

L = -^XibijXj +^2ciXi. (9) 

i,j i 

Now, because the symmetric part of the matrix b{j 
is positively definite (all eigenvalues are positive), 
dL/dt < for sufficiently large \x\. Therefore L ful- 
fils all conditions (2)- (4) of a generalized Lyapunov 
function. 

As an elementary example we consider the Lorenz 
model H 

x = a(y — x); y = —xz + rx — y; z = xy — bz 

(10) 

with r,a,b > 0. Sparrow, see Appendix C in ||, 
proved the boundedness of this model with the aid 
of the following function 

L = rx 2 + ay 2 + a(z — 2r) 2 . (11) 

After the coordinate shift x\ := x, X2 ■= y, X3 := 
z — 2r, we obtain 

L{x) := rx\ + cjx 2 2 + ax\ (12) 

and 

L = -2a(rx\ + x 2 2 + bx\ + 2brx 3 ) (13) 
which is negatively definite for sufficiently large dis- 
tances \j x\ + x\ + x 2 . Thus L(x) as defined in 
(12) fulfils the conditions of a generalized Lyapunov 
function with the implication that a bounded do- 
main exists which attracts all trajectories. 

3 Application to the complex 
Lorenz model 

The complex Lorenz model reads in standard form 
i 

X = - a X + aY 



Y = -aY + rX - XZ (14) 
Z = -bZ + ^(X*Y + XY*) 

where X, Y and Z are complex variables and real, 
respectively. Furthermore, a = 1 — ie, r = r\ + ir 2 
with real parameters e, r%, r 2 , a, and b. In the case 
of modeling a detuned laser, the constants n, r 2 are 
related to the pumping rate and to the detuning, re- 
spectively. Furthermore o = k/j± and b = 711/71. 
where k, 71,711 denote the relaxation constants of 
the cavity, of the polarization, and of the inversion. 
The variable X is proportional to the complex elec- 
tric field amplitude, Y is a linear combination of 
electric field and polarization, which are both com- 
plex, while Z is related to the so-called population 
inversion, for details see e.g. ||. As is well known 
0, this model has nontrivial stationary solutions 
only in the so-called laser case with the parameter 
constraint e = — r 2 . 

It is convenient to introduce real variables Xj, 
with i = 1, 5, by X = xi + ix 2 , Y = X3 + 1x4 
and Z = X5. The real version of (|14|) then reads 

x\ = —ax\ + (TX3 

X\ = —CTX2 + (TX4 

x 3 = r\X\ - x 3 - r 2 x 2 - ex 4 - x\x 5 (15) 

x" 4 = r\X2 - x 4 + r 2 x\ + ex 3 - x 2 x 5 

x 5 = -6x5 + X1X3 + x 2 x 4 . 

In H the following Lyapunov function was pro- 
posed 

L = D 2 {x\ + x\) + x\ + x\ + (x 5 - n - D 2 a) 2 (16) 

which has the Lie derivative 

1 dL t-,2/ 2 , 2\ 2 2 

-— = -aD (x x + x 2 ) - x 3 - x 4 - 

^2(^2^3 - xix A ) - 6x 5 (x 5 - n - D 2 a). (17) 

Here D is an arbitrary parameter at our disposition. 
The latter expression turns out to be negatively def- 
inite for sufficiently large distances \J x 2 + ...x 2 pro- 
vided D obeys the condition 

r 2 /(4D 2 a) < 1. (18) 



3 



After the coordinate shift x' 5 := x§ — r\ — D 2 a, the 
function L{x\,X2, x;j, X4, x' 5 ) := L{x\, ...X5) fulfils 
all requirements (2)-(4) of a generalized Lyapunov 
function. It is thus proved that also the complex 
Lorenz model is bounded for all parameters, with 
and without the laser condition e = — r 2 ||. 

For quantitative results one determines the el- 
lipsoid L = C* which touches the (geometrically 
different) ellipsoid dL/dt = from the outside. 
This amounts to a five-dimensional secular problem 
which in the given case happens to be feasible ana- 
lytically. The attractive domain is then minimized 
with respect to the parameter D with due attention 
paid to the constraint (18). Details can be found in 
@. Numerical evaluations for physically relevant 
parameters give upper bounds for the laser electric 
field which exceed the maximum values reached by 
asymptotic solutions of (15) by factors of between 2 
and 6. In extreme cases of transient evolution, the 
solutions approach within 20% of the upper bounds 
predicted by the Lyapunov method, see Fig. 4 and 
5 in §. 

As a remark, we have examined the more general 
Lyapunov function 

L = D 2 (x\ + x|) + x\ + x\ + 2£ (x 2 x 3 - £1X4) 

v = n + D 2 a + £r 2 ; f < D 2 , (19) 

with the further disposable parameters £ in addition 
to D. L fulfils the properties (2) and (4) as is im- 
mediately seen after the coordinate transformation 

Xl 

X4 
as 



X3, 



xi - iju 2 x 4 , x 2 -> x 2 + i/D 2 x 3 , x 3 
X4, X5 — > xi — v. The Lie-derivative is given 



1 dL 
2~dt 



-(D 2 a + ir 2 )(x 2 +x 2 2 )-x 2 
[r 2 + (cr + l)£](x 2 x 3 - xix 4 ) 



-6x5 + 61^x5. 



(20) 



This derivative turns out to be negatively definite 



for sufficiently large distances y x 2 + ...x 2 , and thus 
obeying (3), provided 



[rg + (a + 
D 2 a + £r 2 



< 4. 



(21) 



The determination of upper bounds of the electric 
field amplitude |AT| 2 is carried out in a similar way 
as in m . The main challenge consists in the task 
of simplifying rather involved analytical expressions 
for different parameter regions. In a physically rel- 
evant parameter domain 



ri < 



4ori(2-b)(2(7-6) 
(a+l) 2 + (2-b)(2a-b)' 



(22) 



with b < 2 and 2a > b, the following upper bound 
is found Q which is minimized with respect to the 
two parameters D and £ 



\X\ 2 < Aan 



(a + 1) 2 + (2-6)(2ct -bY 



(23) 



It is smaller than the upper bound \X\ 2 < \ar\ as 
found previously in || with one disposable param- 
eter only, namely D. 

4 Application to the Navier- 
Stokes equations 

We consider the incompressible Navier-Stokes equa- 
tions (NSE) 



- + (..V). 



-grad(p)+AtAi; + po/ (24) 



in the space domain f2 which, at least in one di- 
mension, has a finite diameter so that it can be 
located between two parallel planes with finite dis- 
tance I. Because of div(v) = 0, the density po is 
constant. We are not concerned here with weak so- 
lutions, which are discussed in H, and assume that 
the solutions are sufficiently smooth. As boundary 
conditions we adopt the no-slip case with v\qq = 0. 

In order to obtain a dynamical system of ordi- 
nary differential equations, we represent the veloc- 
ity field v in terms of an orthonormalized system 
® n G D(Q), with div(<3? n ) = 0, n=l,2,.., where 
D(£l) denotes the space of C°° functions with com- 
pact support in Q. We write 

v(x,t)= J2 c*(*)*n(x) with c n eR (25) 

n=l,2,... 
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and define the Lyapunov function as follows 
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1 dL 2 



|v|| + polM! ||f | 



(29) 



2 dt 
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